Simple Formula for High-Energy Gluon Bremsstrahlung in a 

Finite, Expanding Medium 



Peter Arnold 

Department of Physics, University of Virginia, 
Box 400714, Charlottesville, Virginia 22904, USA 
(Dated: March 5, 2009) 

Abstract 

Previous authors have considered the problem of the medium effects on single gluon 
bremsstrahlung associated with producing a high-energy particle in a finite, time-dependent QCD 
plasma. Working to leading logarithmic order, I show that the result for the bremsstrahlung gluon 
spectrum can be cast into a remarkably simple form in the general case. I similarly analyze the 
process of pair production. Also, I comment on the radius of convergence of the opacity expansion 
in cases where the leading- log approximation holds, showing that the opacity expansion does not 
converge when the thickness of the plasma is greater than roughly the bremsstrahlung formation 
time. Additionally, as a special bonus — available for a limited time only while supplies last! — I 
summarize translations between the notation used by a few of the groups who have worked on this 
and related problems. 



I. INTRODUCTION AND MAIN RESULT 



Roughly a decade ago, Baier, Dokshitzer, Mueller, and Schiff (BDMS) [2] derived a simple 
result for the effect of the medium on the the probability of single gluon bremsstrahlung from 
a high-energy parton produced by some hard process in the background of a uniform, time- 
independent chunk of hot QCD matter (known as a "brick"). Their simple result (based on 
application of a more general formalism) was derived for cases where the number iV co h of 
coherent soft scatterings during gluon bremsstrahlung is large, and they looked for a result 
valid to leading order in (lnA^ coh ) _1 . They found 

d ex 

U du^ 1 ~ Iy3iC ^ = ~ xPs ^ x ) ln l cos (^o£)| , (1-1) 
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where / is the probability of gluon bremsstrahlung from the high-energy particle of energy 
E and species s (quark or gluon), J vac is the corresponding probability had the hard particle 
been produced in vacuum, uj = xE is the energy of the bremsstrahlung gluon, P s _^ g (x) is the 
usual vacuum splitting function, L is the distance the high-energy particle travels through 
the (uniform) medium before abruptly exiting into vacuum, and c^o is a complex number 
with phase exp(— in/ 4) given by 

[(l-x)C A + ^,R 

2x(l-x)E V ; 

Here, C^q is the average squared transverse momentum transfer per unit length that a high- 
energy particle with color representation R picks up through soft, elastic collisions with the 
medium, evaluated at leading-log order, 



q 



where CrT e \ is the collision rate (which is the same at leading order for high-energy quarks 
and gluons, except for an overall factor of the quadratic color Casimir Cr). The leading-log 
approximation arises from the need to cut off the large q± behavior of this integral, which I 
will briefly review later. 

In another paper [3], 1 BDMS showed that they could also find leading-log results for 
non-uniform, time-dependent media, such as an expanding quark-gluon plasma. The result 
was not as simple, however, involving a double integral of a complicated function found 
for the particular case they studied. In this paper, I show that there is a magically simple 
generalization of (1.1) to the general case of non- uniform, time-dependent media. The result 
is 

d ot 

co—(I-I vac ) = -xP s ^ g (x) ln|c(0)|, (1.4) 

du 71 



1 Readers should beware that Ref. [3] investigates a slightly different problem than the one proposed 
here, and gets a correspondingly different answer, for example, for the brick case (1.1). Here, as in Ref. 
[2], I consider radiation from a high-energy parton after it leaves a hard collision that occurs inside the 
medium. Ref. [3], in contrast, purports to study the case where the particle approaches the medium from 
the outside. See the discussion immediately following Eq. (42b) of Ref. [2]. 
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where c(t) satisfies the differential equation 



dt 2 



■u*(t)c(t) 



(1.5) 



with the boundary condition that c(t) approach the constant 1 as t — > oo, and the convention 
that t — is the time of the hard collision that produced the initial high energy particle. 
Here, tVo(t) is (1.2) evaluated at the position of the high-energy particle at time t, and now 
q = q(t,x(t)) is time dependent. The fact that the particle eventually ends up in vacuum 
means that ui^it) — > as t — > oo. 

I will later give the generalization of the result to the case g — > qq of pair production. 

I should note that BDMS's result and my generalization are not complete descriptions 
of the average bremsstrahlung spectrum at leading-log order [4]. For sufficiently small L, 
the average medium effect on bremsstrahlung is instead dominated by atypical events where 
there is a single, larger-than-normal scattering from the medium. I will review this later, 
along with the condition on L [5]. 

The simple form (1.4) is peculiar to three spatial dimensions (i.e. two transverse dimen- 
sions) . I do not know of a generalization that would give a comparably simple result in other 
dimensions. 

In the next section, I review the starting point for the calculation, based on the formalism 
of Zakharov [10] and Baier, Dokshitzer, Mueller, Peigne, and Schiff (BDMPS) [2, 6-8]. I 
organize the notation in a way that's a little friendlier for perturbative calculations in a QCD 
medium with non-static scatterers than the original BDMPS version. (See the discussion in 
the appendix.) Then I review the leading log approximation and the range of validity of the 
BDMS result (1.1). In section III, I derive the basic result (1.4) of this paper. Section IV then 
gives various examples for some cases where the equation (1.5) for c(t) has analytic solutions. 
Section V analyzes the general problem in the limiting cases of a QCD medium that is narrow 
or wide compared to the formation length for gluon bremsstrahlung. Throughout this paper, 
I focus on the case of bremsstrahlung in order to simplify notation, but the formalism applies 
equally well to pair production. In section VI, I give the corresponding results for the case 
of pair production. Finally, in section VII, I comment on implications of BDMS's original 
result (1.1) for the convergence of what is know as the opacity expansion — the expansion of 
the bremsstrahlung probability in powers of the number of elastic scatterings. 

The notational conventions that I use are not exactly the same as those of BDMS or 
Zakharov. The relationship between my notation and various other authors is discussed in 
Appendix A. 

II. STARTING POINT AND ASSUMPTIONS 
A. Notational preliminaries 

Throughout, I will use Cr to denote the quadratic Casimir T^T^ for the color represen- 
tation R associated with some particle, with color generators T^. For a particle of type s, I 
will abbreviate this as C s . For QCD, 



C q = Cp 
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C g = Ca — N c — 3, 



(2.1) 



2N C 
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where N c = 3 is the number of colors, d R will be the dimension of the color representation, 
so that 

rf q = d F = N c = 3, d g = d A = N 2 C - 1 = 8. (2.2) 
tR = Cjid R /dA will be the trace normalization defined by tr(T^T^) = t R 5 ab , with 

t q = *f = \ , t g = t A = N c = 3. (2.3) 

The Dokshitzer-Gribov-Lipatov-Altarelli-Parisi (DGLAP) splitting functions in (1.1) and 
(1.4) are 

P„(.)-ftfi±£l£fl. (2.4) 

Throughout this paper, I will generally place a bar over quantities when I have removed 
an overall factor of Cr from its definition. So I work with q, for example, instead of the 
more standard q. 



B. General Formalism 

Calculations of bremsstrahlung from sufficiently high energy jets must take into account 
the Landau-Pomeranchuk-Migdal (LPM) effect, which arises when the quantum mechanical 
duration (formation time) of the bremsstrahlung process becomes comparable to, or exceeds, 
the mean free time for small-angle elastic collisions. The basic procedure for making such 
calculations was laid out for QED by Migdal in 1956 [1]. The generalization to QCD requires 
accounting for the fact that a bremsstrahlung gluon, unlike a photon, carries (color) charge 
and so can also undergo collisions during the formation time. I will find it convenient to 
start with the particular version of this result derived by Zakharov [9, 10]. This is equivalent 
to the BDMPS formalism of Baier et al. [2, 11], and I will use some of that correspondence 
in how I choose to write Zakharov's result. The general formula is 

d lT v axP s ^Jx) „ f°° , f°° , 

■5"-'-» = w -w R 'l " ti L " h 

VB 1 -yB,{G(B 2 ,t 2 ;B 1 ,t 1 )-G VK (B 2 ,t 2 ;B 1 ,t 1 )}] , (2.6) 

where G(B2, £2! Bi, ti) is the Green's function for a two-dimensional quantum mechanics 
problem with the time-dependent Hamiltonian 

H(t) = 5E(p B ,t)-iT 3 (B,t). (2.7) 

The two terms in H above will be described in a moment. The Green's function 
G(B,t] B\,ti) is a solution to the Schrodinger equation 

idti/>(B,t)=H(t)iP(B,t) (2.8) 

with initial condition 

G(B,t 1 ;B 1 ,t 1 ) = 5 (2 \B-B 1 ). (2.9) 
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The first term in (2.7) describes the energy difference 



(E s P + E gk )-E s p+k * P l±^ + ^5 - '^ + ^ |2 + m - (2.10) 
V SjP T g)fc , s , P +fe 2 ^ -t 2k 2{p + k) V ; 

between (i) a high-energy parton of momentum P = p + k and energy E = P and (ii) 
the same parton with momentum p plus a bremsstrahlung gluon with momentum k. If 
(following Ref. [2]) one defines 

_ pk ± -kp± 

Pb = j , (2.11) 

then this energy difference can be rewritten as 

seta,, t) . + + < {t) 



2x(l-x)P 2(1 -x)P 2xP 2P 
p 2 B + x 2 m 2 (t) + (1 — x) m 2 (t) 
2x{l-x)E ' 



(2.12) 



The notation m(t) accounts for the fact that the effective masses will change as the particle 
transverses a inhomogeneous or time- dependent medium. Qualitatively, the expectation of 
1/SE(pb) is of order the formation time for the bremsstrahlung process in the medium. 2 
The second term in (2.7) is given by 

r 3 (B, t) = \C A T 2 (B, t) + (C s - \C A ) T 2 (xB, t) + \C A f 2 ((1 - x)B, t) , (2.13) 

where f 2 is related to the Fourier transform of dT c i/d 2 q± and defined by 

2 dT cl {t) _j b . q± ,_ 1 f * dT cl {t) _ Jb . q± , 



r 2(M ) ^ J ^ (1 - e-) = - J ^ (1 - e-). (2.14) 

I have not used exactly the same notation as either Zakharov or BDMS, and I summarize 
the differences of notation in Appendix A. On a slightly more substantive matter, both 
implicitly assumed that the rate r e i for soft scattering of the high-energy particle could 
be written as a number density n of static particles in the medium times a cross-section 
(T c i for scattering from such particles. However, their results do not actually depend on 
this assumption. If one simply writes their formulas in terms of the rate r e i rather than 
na, then they apply equally well to the case of scattering from non-static particles, which, 
for example, was analyzed for leading-order calculations in an infinite, time-independent 
thermal medium by Arnold, Moore, and Yaffe (AMY) [12-14] and Jeon and Moore [15]. 
Specifically, the differential rate is 

dT e i s f f d 3 p2 da, 



d 2 q A 



Jd^^J^^fsMil^fs^-q)}, (2.15) 



Here, Cn<j e \ is the soft, elastic scattering rate for a high-energy particle to scatter from a 
single plasma particle of momentum p 2 and species s 2 . q± is the transverse momentum 



2 There is a difference between my use of the phrase "formation time" and Zakharov's [10]. See Appendix 
A. 
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transfer to the high-energy particle from this single scattering. f{p2) is the phase space 
density of plasma particles per degree of freedom, which in thermal equilibrium is the Bose 
or Fermi distribution for the plasma particle. v S2 is the number of spin, color, and flavor 
degrees of freedom for species s 2 , which would be 2d a = 16 for gluons and Ad F N { = 12N { for 
the sum of quarks and anti-quarks, where JVf is the number of quark flavors. The factor of / 
gives the density of plasma particles, while the factor of 1 ± / is a final-state Bose or Fermi 
factor. Final state factors for the high-energy particle (as opposed to the plasma particle 
it is scattering from) may be included at the end of the LPM calculation, if desired (see, 
for example, the 1—^2 splitting terms in the effective kinetic theories of Refs. [12, 13, 16]), 
but in the present context I assume that the high-energy particle is an isolated particle of 
energy much higher than the plasma particles, so that its final state factor can be ignored. 

In terms of specifics, perturbative calculations for a QCD plasma in local equilibrium 
give the simple formulas 3 



dT, 



el 



g 2 Tml 



d 2 q± (2tt) 



x < 



q 2 ± (ql + ml) 



q± « T, 
q± » T, 



(2.16) 



in the limits of q± small or large compared to the temperature T. Here mo is the Debye 
mass, 

ml= (t A + N { t F y 3 g 2 T 2 = (l + |iV f ) <? 2 T 2 , (2.17) 
and M is the weighted number density 

* s £ / *•><*) = §§ («* + r = §§ ( x + &<) T "' < 2 - 18 > 

where ((z) is the Riemann zeta function. 

The formalism reviewed above assumes that the characteristics of the medium do not 
change significantly over a Debye screening length. It is not restricted to equilibrium situa- 
tions, but I will assume that the differential elastic cross-section is isotropic in the transverse 
plane. The formalism also assumes that the final bremsstrahlung gluon and accompanying 
particle are energetic enough that transverse momentum transfers from the medium will be 
small compared to their momenta. 

More generally, all calculations based on variations of Migdal's procedure require that the 
mean free path for soft, elastic collisions be large compared to the screening length. 4 This 
assumption holds for a thermal plasma in the weak coupling limit, where the mean free path 
is order l/g 2 T and the screening length is order 1/gT. (See, for example, the discussion in 
Ref. [12].) 



3 The simple form of the q± <C T formula comes from Ref. [17]. This is the formula used by AMY [14] in 
studying the LPM effect in hydrodynamic transport coefficients, where the relevant particle energies are 
E ~ T. 

4 More precisely it is the mean free path for the subset of soft, elastic collisions which contribute to the 
result at the desired accuracy. In the thermal case, for example, ultra-soft magnetic interactions with 
Q± ~ 9 2 T do not affect results at leading order in coupling. 
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If running of the coupling constant a is included in the analysis, then dT c \/d 2 q± should 
plausibly be evaluated with g 2 (q±). 5 In Ref. [19], it is argued that the overall factor of 
a associated with the coupling of the bremsstrahlung gluon [here the explicit a in (1.4) or 
(2.6)], should plausibly be evaluated as a(Q±), where Q± is the typical transverse momentum 
transfer over the formation time and is discussed below. This last prescription is in the spirit 
of earlier suggestions by BDMPS [8]. 6 

C. Leading Log (Harmonic Oscillator) Approximation 

Consider a medium that is thick enough that the total number of soft scatterings with 
individual momentum transfers > m D , as the particle traverses the medium, is large. 7 
In the high energy limit, the number iV co h of such scatterings in a formation time also 
becomes large. As noted long ago by Migdal [1], the calculation of the LPM effect simplifies 
significantly if one works to leading order in (mA^ch) -1 . In the context of QCD, such 
leading-log calculations were explored by Baier et al. using their BDMPS formalism and 
what is known as the harmonic oscillator approximation. Following BDMS [3], I will focus 
on leading-log calculations in this paper as well. 

The large N co ^ limit corresponds to the case where the total transverse momentum trans- 
fer Qj_ to a high-energy particle during the formation time is large compared to the screening 
mass m-Q. One consequence of large transverse momentum is that we can ignore the effective 
particle masses m s and m g in (2.12). Another consequence is that large Q± corresponds in 
Fourier space to small B. Naively, Eq. (2.14) for F2 can then be replaced by its small b limit, 
which is formally 



This is known as the harmonic oscillator approximation because of the form of (2.19). The 
problem is that the above integral is logarithmically divergent because of the large q± be- 
havior of (2.16). For a leading log analysis of typical events, it should be cut off at order of 
the typical total momentum transfer Qj_ in a formation time. Parametrically, recalling the 



5 See, for example, Rcfs. [18, 19]. In order to avoid an unphysical infrared divergence of the calculation 
when the definition of g 2 blows up at Aqcd , one should appropriately cut off the running in the infrared. 
One possibility would be to use g 2 {\Jq]_ + ti d ). 

6 Specifically, after Eq. (3.12) of Ref. [8], they suggest taking a s (k) with k oc L 1/>2 for the calculation of 
average bremsstrahlung energy loss in a thin QCD medium. For that problem, the energy loss is dominated 
by gluons whose formation time is of order the length L of the medium. In that case, Q±_ oc L 1 / 2 as in 



7 This statement contains the restriction q± > mD because, in the weak coupling limit, the most common 
scatterings, by a parametric factor of ln(a: _1 ), have momentum g 2 T < q± <C and are mediated by the 
exchange of low-frequency magnetic gluons. These low-frequency magnetic gluons are not Debye screened, 
and their contribution is cut off only by non-perturbative effects. However, these ultra-low momentum 
scatterings do not contribute at leading order to (1.3) and (2.14) [because of the factor of q\ in (1.3)] and 
so do not have an effect on bremsstrahlung at leading order in coupling. 




(2.19) 



(2.20) below. 
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definition of q, 



.(C Rt qL) 1/2 , L<L CX) 



where L is the characteristic thickness of the medium and L cr is the infinite-medium forma- 
tion time 8 

/ E \ 1/2 

Above, Ei is the energy E, xE, or (1 — x)E of a particular parton in the splitting process, 
and one should use whichever parton gives the smallest Q±. For small x, that will be the 
bremsstrahlung gluon, giving L CT ~ (uj / C^q) 1 / 2 . 

Using (1.3) and (2.16), the leading-log value of q for a weakly-coupled thermal QCD 
plasma is then 9 

f~aTmjUn(%) (2.22a) 
\ron/ 



if Q± < T and 



ml J \ T 2 



)2 



aTm'l In ( ^] + 4na 2 Af\n ( ) (2.22b) 



otherwise. For 3-flavor QCD, aTml and ArKa 2 M differ by only about 15%, and so one 
could combine the logarithms of (2.22) into either aTml \n(Q 2 ± / ml) or An^Afln^Q^/ml) 
without much error. 

If Q± is so large that a(Q±) is significantly different from a(m D ), then one should include 
1-loop running of the coupling when integrating (2.16). The result can be put into the form 10 

q-~a(T)Tml(m D ) In (J^j + An a(Q ± ) a(T) Af\n (Q^j , (2.23) 

where m^mjy) indicates the Debye mass (2.17) evaluated with running coupling g 2 (mv). 

Note that the leading-log formula (2.22) for q depends logarithmically on Q±, which in 
turn depends on q. One could determine q self-consistently, but it should be kept in mind 
that a precise value of q inside the logarithm is not called for because we are only pursuing 
a leading- log result. For an example of how things work out at next-to- leading logarithmic 
order, see the infinite- medium calculation of Ref. [19]. 

In any case, in the leading- log approximation (2.19), the 2-dimensional Hamiltonian of 
(2.7) becomes a 2-dimensional harmonic oscillator problem with time-dependent frequency: 

H(t)^^L + lMu 2 (t)B 2 , (2.24) 



8 See, for example, the discussion in Sec. 3 of Ref. [7]. 

9 For (2.22a), see also Eq. (13) of Ref. [20] and the relation to Ref. [21] discussed after Eq. (61) of Ref. [20]. 
See, for example, the discussion in Sec. VI of Ref. [19]. Though the form of (2.23) is convenient, it can 

be misleading. In the limit that Q± is so large that a(Q±) -C a(T), the answer does not actually depend 
on q± of order Q± — it is instead dominated by those q± for which a(q±) is of order a(T) [19]. Also, the 
simple formula (2.23) is only valid if there are no vacuum mass thresholds between win and Q±. 
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with 



M = x(l-x)E, 



ul(t) = -i 



[i- X )c A +x 2 c s ] m 

2x(l-x)E 



(2.25) 
(2.26) 



Note that u 2 is imaginary. Its inverse magnitude l/|c^o| is of order the infinite-medium 
formation time L cr of (2.21). 

The harmonic oscillator approximation breaks down for sufficiently small L, even when 
logarithms are large. Using the typical total momentum transfer (2.20) as an upper cut-off to 
determine the integral in (2.19) ignores the possibility of bremsstrahlung from rare, atypical 
scatterings with lager q±, which turn out to be important for sufficiently small L. There 
has been some confusion about the resulting range of validity of the harmonic oscillator 
approximation used by BDMS for a leading-log analysis of the spectrum. Zakharov [4] 
suggested that the harmonic oscillator approximation outlined in this section is only valid 
when L ^> L cr , which is equivalent to \ujqL\ ^> 1. In Ref. [5], however, I argue that the 
validity extends to 

(2.27) 



[HQl/ml)] 1 ' 2 ' 



This includes the interesting region L ~ L CT (equivalently \uqL\ ~ 1) in a leading-log analysis, 
which treats the logarithm as large. 



III. DERIVATION 



A. A double integral 



If G is the Green's function, then the two components of the vector function 
W s 1 G(B,t; Bi,ti) will also satisfy the Schrodinger equation (2.8) but with initial condition 



V Bl G(B, h; B 1 ,t 1 ) = V B / 2) (B - B ± ). 



(3.1) 



In (2.6), we are interested in the particular case Bi = 0, which then corresponds to the 
initial condition 



V ' b 1 G{B, ti, B 1 , ti 



Bi=0 



-V B 5 (2 \B). 



The desired solution in the leading log approximation (2.24) is 



V Bl G(B,t; B 1 ,t 1 ) 



M 2 



Bi=o 2nS 2 (t;t 1 ) 
where S(t; t±) satisfies the differential equation 

d 2 S = -u 2 (t)S 

with boundary conditions 

S(t 1 ;h)=0, dtSfah) 



B exp 



iMd t S(t; t 



— B : 



t=ti 



(3.2) 



(3.3) 



(3.4) 



(3.5) 
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One may check this by (i) plugging it into the Schrodinger equation and noting that it is 
a solution, and (ii) checking the initial condition by solving for t infinitesimally close to ti, 
where S(t; ti) — > t — t\ and (3.3) becomes 

M 2 „ ( %M 
B exp 



'VbjG^B, t; B 1 ,ti) 



Bi=0 



2Tt(t-t 1 ) 2 

M 

-V B ; 



exp 



2(t-h) 
M 



B' 



'2ni(t -h) 

which is a representation of —'Vb^(B) for infinitesimal t — ti 
Substituting (3.3) into (2.6) gives 



2i(t-h) 



d 



a 



— x P s ^Jx) ReX, 

7T 



where 



1= - 



dti 



dt? 



(3.6) 



(3.7) 



(3.8) 



-S 2 (t 2 ;ti) (t 2 -n, 

This gives an answer in terms of a double integral involving the function S(t 2 ] ^i)- But both 
integrals can be done explicitly, even for the case of arbitrary u>o(t). 



B. The t 2 integration 



Now consider the other solution to the 2nd-order differential Eq. (3.4), which I will call 
C(t;ti) and take to have boundary conditions 



C(t 1 ;t 1 ) = l, dtC{t-M 



0. 



t=ti 



(3.9) 



If ui 2 {t) were a constant, then the two solutions would be S = cUq 1 sin(uj (t — ii)) and 
C = cos[cuo(t — ti)), which is the motivation for the labels S and C. 
The form of the differential equation implies that the Wronskian 

W = Cd t S -Sd t C (3.10) 
is independent of time and so always equal to its value at t — ti. 

Cd t S- Sd t C = 1. (3.11) 
Dividing both sides by S 2 then gives 



dti^)=- 2 . 



(3.12) 



We can use this to do the t 2 integral in (3.8). Rewrite the time integrals in (3.8) to have 
upper limit t, taking the limit t — > oo at the end. Then rewrite the t 2 integral as 



lim / dt 2 



S 2 (t 2 ;h) (t 2 -h) 



So 



X = lim 



dti 



C(t;ti) 
~S(t;h) 

~S(t;h) 

Cjt-h) 

S(t;t x ) t-h 



1 

h lim 

* - h e^O 
1 



C(ti + e;ti) 1 
S(ti + e;ti) e 



1 



(3.13) 
(3.14) 
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C. The t\ integration 

Now note that 

C(t;t 1 ) = -d tl S(t;t 1 ). (3.15) 

This follows because (i) —d tl S will satisfy the same equation (3.4) that S does, and (ii) the 
boundary conditions work out correctly. The boundary conditions (3.9) can be confirmed 
from the small t — t\ expansion of S(t; ti), which is 

S(t; h) = (t- tl )-^ u>l{h) (t - ti) 3 + 0[(t- h) 5 ], (3.16) 

so that 

C(t; t,) = 1 - 1 ulih) (t - t,f + 0[(t- t x ) 4 ]. (3.17) 

I shall not need it, but the corresponding derivative of C is 

uj 2 (t 1 )S(t;t 1 ) = d tl C(t;t 1 ). (3.18) 

Note that the relations (3.15) and (3.18) involve t\ derivatives — I will discuss the case of t 
derivatives later. 

Now substitute (3.15) into (3.14): 

X = — lim In 

Combining with (3.7), 

u-^(I-I yac ) = -xP 9 ^ e (x) limln|W;0)|. (3.20) 

did TT t->oo 



g(Mi) 
t - h 



Ll=Z 



ti =n 



lim In 



S(t;0) 



\im\n[d t S(t;0)}. (3.19) 



D. Final simplification 

The result (3.20) is perfectly adequate, but it is amusing to put it in a final form that is 
even more closely analogous to the result (1.1) for the brick problem. 

Note that any solution to a linear differential equation can be written as a superposition 
of others. So S(t] ti) and C{t; ti) can be expressed as superpositions of S(t; t ) and C(t; t ) 
for any to- Specifically, 



S(t; h) = C{h- 1 ) S(t; t ) - S(t i; t ) C(t; t ), 

C(t; t ± ) = -d tl C(tu to) S(t; t ) + dt.Sih; t ) C(t; t ). 



(3.21) 
(3.22) 



To verify these formulas, one just needs to check the boundary conditions. The conditions 
S(ti,ti) — and d t C(t; ti)\ t=tl =0 are easy. The other two, d t S(t; h)^^ = landC(ti;ti) = 
1, follow from the time independence of the Wronskian, 



C(f ; t ) d t S(t; t ) - S(t; t ) d t C(t; t ) 



J t=ti 



C(t;t )d t S(t;to) - S(t;t )d t C(t; t ) 



I. 



-I t=t 



(3.23) 
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From (3.21), we see that S(t;ti) is anti-symmetric in its arguments: 

S(t 2 ;t 1 ) = -S(t 1 ;t 2 ). (3.24) 

We can then combine this with (3.15) for the t 1 derivative of S to get a formula for the t 
derivative: 

d t S(t; t ± ) = -dtSih; t) = C{h; t). (3.25) 

[Eq. (3.22) does not allow us to deduce any comparable symmetry property of C] We can 
now use (3.25) to rewrite (3.20) in the form 

d ol 

to— (/-/vac) = -xP s ^ g (x) ln|C(0;oo)|, (3.26) 

CtUJ 7T 

which is Eq. (1.4) of the introduction. 



IV. EXAMPLES 



One can of course solve the differential equation (1.5) numerically for any desired time- 
dependence of q(t) along the path of the particle. In this section, I give a few examples that 
have analytic solutions. 



A. The brick problem 

Consider the case where the particle travels distance L through a uniform medium and 
then emerges into vacuum. So 

"*« - {f \ : t ^ 

The solution c(t) to (1.5) is then 

fcosML-f)), t<L; 

\l, t>L. K J 

Eq. (1.4) then reproduces the result (1.1) of BDMS [2]. 

Using the fact that cu is proportional to (— i) 1 / 2 , one can alternatively write the result 
solely in terms of real quantities using the identity 

In | cos(e"^ /4 a;)| = \ In \\ cosh( v / 2a;) + \ cos(v / 2x)j . (x real) (4.3) 

The large L behavior is 

In | cos(cj L) | ~ - In 2, (4.4) 
v2 

up to exponentially small corrections. (But you shouldn't take seriously the In 2 term because 

remember that I've only treated u itself up to leading-log order.) In this limit, one can 
write 

u±{I-I^)^J^ L, (4.5) 
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with 

For fixed x, the small L behavior is 11 

ln|cos(cu L)|~^(|cu |L) 4 . (4.7) 

Small L in this context means \ouo\L <C 1, equivalent to L L cr . But keep in mind that 
the harmonic oscillator approximation breaks down for calculations of the spectrum when 
L < LJiHQl/ml)}^ [5]. 

B. Exponential Profile 

Consider an exponential profile 

u; 2 (t) = u; 2 (0) e-^. (4.8) 

The solution is 

c(f) = J (2o;o(0)Le- t / 2L ), (4.9) 

giving 

4(/-U = -^ s » ln|Jo(2w (0)L)|. (4.10) 

did IT 



C. Power Law Relaxation 



Motivated by modeling Bjorken expansion, BDMS [3] considered the case where q falls 
like a power of time and then suddenly vanishes (the particle emerges into vacuum) at time 
L. So 

W S(t) = H(*b)(W B , t <t<t + L; 
0W \0, t + L<t; K J 

where a is some power and I've now labeled the time of the initial hard process as to rather 
than zero. The solution to (1.5) is then c(t) = 1 for t > t + L (the vacuum solution) and 

/.v f Z Y J "( z ) Y v-i{zl) ~ Y v (z) J v -i(z L ) 
c{t) = — 



zlJ Jv{zl) Yv-i(zl) - Y u (z L ) J v -i(z L ) 

(j-^j [Uz)Y u ^{z L )-Y u {z)J u ^{z L )\ (t<t + L) (4.12) 



irz L I z • ' 



Readers familiar with the fact that the medium-induced contribution to energy loss is proportional to qL 2 
for small L [7] may wonder how the L 4 behavior in the spectrum (4.7) is consistent. In (4.7), the limit is 
that L is small compared to the formation time, which is of order (xE/q) 1 / 2 = {uj/q) 1 / 2 for x not close to 
1. In contrast, the small L formula for energy loss assumes L <S^. (E/q) 1 ' 2 . In the latter limit, when the 
energy loss is determined by integrating udl/du over o>, the integral is dominated by o/s for which the 
formation time is of order L (u <~ qL 2 ), where the small L assumption of (4.7) has just started to fail. 
Using (4.7) merely as a parametric estimate then yields AE ~ au\uJo\ 4 L 4 ~ auj(q/uj) 2 L 4 ~ aqL 2 . 
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where 12 



The final result is then 



1 

2-a' 

z = z(t) = 2vu (t ) t ^j- 

z = z(t Q ), 
z L = z(t + L). 



l/2v 



UJ^~{I - / vac ) = -XP S ^ S (X) In 

auj it 



to + L 



1/2 Jv(zq) Yv-\{zl) ~ Y V {Z G ) J„-i(z L ) 



Jv(zl) Yv-i(z l ) - Y u (z L ) J v -x(z L ) 



(4.13) 

(4.14) 

(4.15) 
(4.16) 



(4.17) 



D. sech 2 Profile 

As a final analytic example, consider a hard particle starting at t — t with profile 



The solution is 



uj (t) = VI sech 



- ^ac) = -xP s ^ g (x) In F(a+,a„; 1; 
auj 7i 



e 2t /L _|_ 



t) 



where F is the hypergeometric function and 



a± = Uyi + (2fii:) 2 . 



(4.18) 



(4.19) 



(4.20) 



V. GENERAL SOLUTION: LIMITING CASES 

I now turn to the behavior of the general solution (1.4) for the limits of small or large 
width of the medium for fixed x. 



A. Small width 

In this case, we can solve the differential equation 

c(f) = -ul{t) c(t) 



(5.1) 



12 My z differs by a factor of i from that of Rcf. [3], which is why they have modified Bessel functions K 
and I instead of J and Y. Also, an equivalent way of writing (4.12) is to replace Y v _\ and Y v by Ji_„ 
and J_„. If comparing to Ref. [3], keep in mind that they solve a slightly different problem, as explained 
in footnote 1. 



14 



by perturbing around the vacuum solution c vac (t) = 1. The solution is 

c{t) = l + Cl (t) + c 2 (t) + O(cu 3 ), (5.2) 

where 

POO 

Cl (t) = - / dt'(t' - f)wg(t'), (5.3) 



/oo 
dt\t' -t)ul{t') Cl {t'). (5.4) 

Now recall that oJq is proportional to — i so that ci is imaginary and c 2 is real. Then 
ln|c(0)| ~|ln[(l+c 2 (0)) 2 +| Cl (0)|' 



~ 2 



2 



|| Cl (0)| 2 + c 2 (0) 



OO \ ^ POO POO 

dt' t' \ul{t')\ - / dt' t'\uol{t)\ / dt" (t" - t') \uol{t")\. (5.5) 
J Jo Jt> 



This is the general form of the small-width answer, of which (4.7) is a specific case. 
B. Large width 

Now consider the case where cJqCO i s a very slowly varying function of t. Then we can 
make an adiabatic approximation, and the most important feature of the solution for s(t) 
will be a "phase factor" that is approximately 13 

c(t) ~ exp J dt' uj (t')^j = exp jf" dt' \uj (t')^j exp dt' \co (t')\^ ( 5 - 6 ) 

Neglecting pref actors (whose effect is parametrically smaller than the exponent), 



1 



5C 



ln|c(0)|~-^ / eft MO | . (5.7) 



o 

Comparing to (4.6), this gives 

(l) ±{I-I yaa )=j^dtu^{t), (5.8) 

as you would expect: In the limit of very thick, slowly varying media, you just treat the 
problem as locally uniform, use the result for the bremsstrahlung rate in an infinite, uniform 
medium, and integrate. 



13 Because ujq is proportional to exp(— mt/4), the other solution exp (— i J°° dt'u>o(t')) is, in the large width 
limit, exponentially small at t = 0, and so its contribution to c(0) can be neglected. 
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There is a technical subtlety if one goes to next order in the adiabatic expansion and 
looks at the prefactor. The assumption of the adiabatic expansion is that \uo\ <C To 
first order in the prefactor, the solution for C{ti,t 2 ) is 



C(t;t 2 ) 



u (t) 



1/2 



exp 



dt' u (t') 



(5.9) 



Then 



In |c(0)| 



1 



dt \wo(t)\ 



In 



^o(oo) 



^o(O) 



(5.10) 



Since cj (oo) = 0, this answer suffers from a logarithmic divergence. 

The problem is that the adiabatic assumption |tL>o| ^ l^ol mus t break down at sufficiently 
late times. As an example, consider the exponential distribution of Sec. IV B. The adiabatic 
assumption first breaks down when |u;o| drops to \uo\ ~ 1/L. If we use this value of ujq to cut 
off the logarithm in (5.10), then we find a correction to the bulk result of size — | ln(|u;o(0)|L). 
And in fact, exactly such a correction appears in the large L expansion of the exact result 
in (4.10), which gives 



In | Jo(2cu (0)L) | = V2\u (0)\L - \ ln(|wo(0)|L) + 0(1). 



(5.11) 



VI. PAIR PRODUCTION g -> qq 

Previous results are easily modified for the case of pair production g — > qq. First, one 
uses the appropriate DGLAP vacuum splitting function, so the overall result (1.4) becomes 

d (y 

x—(I - /vac) = -xP g » In |c(0)| , (6.1) 

dX 7T 

with 

P g ^ q (:r) = N { t F [x 2 + (1 - xf] (6.2) 

if one sums over all quark flavors. Here x is the momentum fraction of the quark. One 
must also change the factors in the definition (1.2) of Wq, as I shall discuss. The only other 
change necessary is to appropriately change the group factors in Eq. (2.13) for T 3 to reflect 
the different arrangement of color representations in the splitting process from F — > AF to 
A -> FF. The generalization of (2.13) is 14 

r 3 (B,t)= i(O i?2 +C R3 -C i?1 )f 2 (a; lJ B,t) 
mC R3 + C Rl -C R2 )T 2 (x 2 B,t) 

mC Rl +C R2 -C R:i )T 2 (x 3 B,t) (6.3) 
for ai?i-> -R2-R3 splitting process with corresponding momentum fractions 

x± — 1, x 2 = x, x 3 = 1 — x. (6.4) 
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For this form, see the discussion surrounding Eq. (6.11) and footnote 24 of Ref. [12]. 
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For s — > gs processes, this gives (2.13). For g — > qq, the color factors of (2.13) (or equiva- 
lently the momentum fractions) are permuted to 

T 3 (B, t) = (C F - \C A ) f 2 (B, f) + \C A T 2 (xB, t) + \C A f 2 ((1 - x)B, t) . (6.5) 

The resulting value of uj 2 replacing (2.26) is then 

2 .[C F -x(l-x)C A ]q 

Ofi = —I ; r . 6.6 

2x(l-x)E v ; 

VII. CONVERGENCE OF THE OPACITY EXPANSION 

The opacity expansion investigated by Wiedemann [22] and Gyulassy, Levai, and Vitev 
(GLV) [23] involves analyzing bremsstrahlung in the QCD medium by expanding order by 
order in the number of elastic scatterings. It is interesting to ask what happens if such 
an expansion is made in a case where the leading-log calculation of BDMS is valid. An 
expansion in powers of elastic collisions is equivalent to an expansion in powers of T 3 (2.13), 
which in leading-log approximation is equivalent to an expansion in powers of u 2 (2.26). 

Now consider BDMS's result (1.1) for the brick problem, and rewrite it in the form 

oo^-(I-I vac ) = —xP s ^ g (x) ln[cos(e^ / V/ 2 )cos(e-^ /4 2 1/2 )l, (7.1) 



duj ' 2% 



where 



z = \ujI\L 2 . (7.2) 
The opacity expansion of this result is its Taylor series in z, proportional to 



In 



cos(eWV/ 2 ) cos(e-/V/ 2 ) = \z 2 - ^ + - • • • . (7.3) 



Mathematically, the expression (7.1) is an analytic function of z, and therefore its radius 
of convergence is given by the distance to the nearest singularity in the complex z plane. 
The nearest singularities are the branch points of the logarithm where either of the cosines 
vanish, at z — ±i(n/2) 2 . In this example, the opacity expansion therefore only convergences 
for \z\ < (vr/2) 2 , which corresponds to 

L < y^- . (brick) (7.4) 
l^o I 

Recall that, qualitatively, l/|<^o| is of order the formation time. The conclusion is that the 
opacity expansion does not converge when the medium is thicker than roughly the formation 
time. 

In Fig. 1, I show the function (7.3) vs its expansion to nth order in the opacity expansion 
for several n. One can see the failure of convergence beyond z = (tt/2) 2 . 

One of the uses of the opacity expansion has been as a hook to derive general results 
by summing up the expansion to all orders, arriving at formalism related to BDMPS and 
Zakharov (for example, as in Ref. [22]). In this case, the lack of convergence of the Taylor 
series for large L does not matter. 

Readers may wonder at the juxtaposition of the opacity expansion and the leading loga- 
rithm approximation. In the large iV co h limit of the leading logarithm approximation, z — 1 
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FIG. 1: The function of (7.3) [solid line] vs z = |u;o| 2 L 2 compared to its Taylor series expansion to 
nth order for selected values of n [dashed lines] . 



in Fig. 1 represents a very large number of elastic scatterings. But the answer is nonetheless 
reproduced well by the n=4 curve, which only includes up to four scatterings. How can 
this be? The reason is that the LPM effect causes even a large number of scatterings to 
behave like a single scattering if they occur within a distance small compared to the forma- 
tion time. For this reason, it is possible for just four scatterings, spread out across L, to 
reproduce the same total bremsstrahlung rate as a large number of scatterings, in leading 
log approximation. 

I should clarify that the expansion discussed here depends on first making the leading log 
approximation, treating q as a constant, and only then making the opacity expansion. So, for 
instance, I have ignored the fact that the upper limit Q± of the logarithm in (2.22) depends 
on the number n of collisions. In particular, readers familiar with the opacity expansion 
may wonder at the absence of a leading n — 1 term in the expansion (7.3), proportional to 
z. This is a special consequence of the leading-log approximation [4, 5]. 15 

One might wonder whether the lack of convergence is an artifact of the brick problem, 
where oo^t) is not an analytic function of time. However, one can draw the same conclusion 
from the exponential profile (4.10). In this case, the singularity occurs at the first zero 
of the Bessel function, when its argument is 2.40482 • • • . The corresponding condition for 
convergence of the opacity expansion in this case is 



L < 



2.40482 



2\u 



(exponential) 



(7.5) 



Specifically, consider Eq. (6.7) of Ref. [22], using definitions (3.39), (5.6-9) and (5.11) of that reference. 
The leading-log approximation is a(p) oc p 2 , which corresponds to £(<Z_i_) oc V 2 ^ 2 -* (q±). If one uses this 
form of E and integrates Eq. (6.7) of Ref. [22] over all bremsstrahlung gluon transverse momenta k± 
(making the k± <C k approximation by integrating all the way up to k± = oo), one finds a zero result. 
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These results have been derived in the leading-log approximation. In situations where 
corrections to the leading-log approximation are small, one expects similar conclusions since 
small perturbations will not remove the presence of singularities. The non-convergence of the 
opacity expansion might possibly be related to the observed poor convergence in numerical 
results at small x by Wicks, shown in Appendix B of Ref. [27], since Lj ' L CT oc x^ 1 ^ 2 in the 
small x limit. 
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APPENDIX A: RELATION OF NOTATION TO OTHER AUTHORS 
1. Zakharov 

The equations I give in Sec. II B are organized slightly differently than by Zakharov 
[10]. I will compare my conventions specifically to Ref. [10]. There, Zakharov ignores the 
fact that effective particle masses in a non-uniform medium will depend on position and 
time. If one treats them as constant, then their contribution to the Hamiltonian defined 
by (2.7), (2.12) and (2.13) is an additive constant, and their sole effect is to contribute a 
simple phase exp[— i(constant) At] in the Green function, which Zakharov explicitly factors 
out. Specifically, the relationship between my Hamiltonian and Green function and those of 
Zakharov (Z) [10] is 



ff(my) — Jj(Z) _|_ 



r(Z) ' 
L f 



G (m y\B 2 ,t 2 ;B 1 ,t 1 )=exp 



i(t 2 



r(Z) 



K^(B 2 ,t 2 \B 1 ,t 1 ), 



where 



2x(l - x)E 
x 2 m 2 + (1 — x) m 



(Al) 
(A2) 

(A3) 



is what Zakharov calls the formation time. Zakharov chooses to incorporate 1/Lf into his 
Hamiltonian in the later work of Ref. [11]. 

There is a difference between his and my use of the phrase "formation time." Zakharov 
uses it to mean the formation time in vacuum in the case of exactly collinear bremsstrahlung, 
which is given by the inverse of (2.12) with ps set to zero. I use it to mean the formation 
time of typical bremsstrahlung in the medium, consistently accounting for the LPM effect, 
which is the inverse of (2.12) including the expectation of p 2 B . 

The gradients in (2.6) correspond (up to factors of +i and —i) to the operators p in 
Zakharov's definition of <?(£i, £2, 
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Finally, the way I have written Zakharov's three-parton and dipole cross-sections 03 
and 02 can be taken from BDMS's discussion of the equivalence of BDMPS and Zakharov 
formalisms in Ref. [2], as I shall discuss below. 



2. BDMS 



Throughout, where BDMS [2] or the earlier works of BDMPS [6-8] expresses rates in 
terms of density p times a cross-section a, I instead write a rate T. This allows one to more 
easily apply the formulas to calculations that account for the dynamical nature of screening 
in the plasma. 

BDMS and BDMPS characterize the differential elastic cross-section in terms of a nor- 
malized quantity 



2 1 da e \ 



V(Q) = (A4) 
<7 e i d z Q 

where they define Q = q±/m D . Translating to the language of rates, one may equivalently 
write 

nQ 2 ) = ^% (as) 

where r c i is written A -1 = pa in the BDMPS formalism and A is the mean free path for elastic 
collisions. This expression is problematical for full, leading-order perturbative calculations, 
however, because the total elastic scattering rate r e i for a high-energy parton traveling 
through a QCD plasma has a logarithmic infrared divergence in perturbation theory, as can 
be seen by integrating (2.16) over d 2 q±. [The divergence does not appear in the discussions 
of BDMS and BDMPS because, when they specialize to the case of Coulomb scattering, 
they model V(Q 2 ) as proportional to \/{q\ + mf^) 2 rather than the actual low-momentum 
perturbative behavior of (2.16).] The divergence arises from the exchange of low-frequency 
magnetic gluons, which are not screened, and is cut off only by the non-perturbative physics 
of magnetic confinement in hot QCD at a momentum scale q± ~ g 2 T. Formally, it is not 
clear whether there is any rigorous, convention independent, non-perturbative definition of 
the total rate A -1 = r e i, and so it is best to avoid the quantity altogether. Fortunately, this 
is just an issue of normalization convention. The various quantities in the BDMPS formulas 
for the bremsstrahlung rate appear in combinations where the factors of A cancel, and I have 
chosen to avoid them in the formulas of Sec. II B. 

A notational translation table is provided in Table I. The reasons for the complex con- 
jugation that appears in some entries of the BDMS column is that BDMS pick conventions 
where their analog of the Schrodinger equation (2.8) corresponds to a problem with nega- 
tive mass M. One can convert to a Schrodinger equation with a positive mass (Zakharov's 
convention, which I adopt) by taking the complex conjugate of the equation, which takes 
V> -> ip*, M -> —M, and u -> u%. 

In BDMS [2], the quark and gluon masses are ignored. This is parametrically valid when 
Q± ^> m q and m g , which for a thick medium (L > L CT ) corresponds to the high-energy limit 
E ^> rrift/ q. In perturbation theory, where m q ~ m g ~ gT, this condition is parametrically 
E 3> T. However, in applications of the LPM effect where E ~ T is of interest (such 
as leading-order calculations of viscosity and other transport coefficients [14]), one should 
include the mass terms. 

Finally, there is an overall minus sign difference between my (2.6) and the comparable 
Eq. (59) of Ref. [2]. One quick way to resolve minus sign issues is to check that the final 
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this paper 


BDMS [2] 


Zakharov [10] 


AMY [12] 










<?_L 






<?± 


fcl 


1 per 












dr c i 

<Pq ± 


V(Q 2 ) 
/i 2 AC F 




J2^J 2n\ A \Q)\A \Q)\ ) q o= q z 


PB 




P 


h/p> 


B 




P 




<f 


C F A 
u>£N c 


|nC 2 (0) 






2C F \ 






t 


2C f Xt 






M 


kN c _ 2ACf// 2 
2AC f ~~ Ar c 


/i(rr) 




G(B2, i 2 ; Bi,ti) 


f x 2 G*(B 2 ,z 2 ;B 1 ,z 1 ) 


e -««2-«l)/i/ 

xX(p 2 ,6|Pl,?l) 




f 2 (B) 


1-V(B) 
AC F 


|na 2 (p) 





TABLE I: Translation between notation of this paper and various authors. The entries in the q line 
are logarithmically divergent in the ultraviolet and should be understood as appropriately cut off 
for a leadingdog approximation, as discussed in the text, or equivalently evaluated at some small 
effective value of B (BDMS) or p (Zakharov) of order 1/Q±. 

answer for the effect of the medium is positive in the limit of a very thick medium, as in 
(5.8). 16 



AMY 



Next, I wish to make contact with the notation used in my previous work with Moore 
and Yaffe [12-14]. That analysis was for the case of an infinite, uniform, time-independent 
medium. Following Migdal [1], one can treat this case by starting with the non- vacuum part 
of (2.6), changing integration variables from t 2 to the time difference At = t 2 — ti, and then 
using time invariance to note that the Green function depends only on At. The t\ integral 
then just gives a factor of the total time, and the resulting equation for the bremsstrahlung 
rate is 



Now define 



b 1= b 2 =o 



f{B,t) = 2i V Bl G{B,t;B 1 ,0) 



J Bi=0 



(A6) 



(A7) 



16 There appears to be a lost minus sign in the transition from Eqs. (31) and (33) to (51) of Ref. [2], which 
then propagates to their (59). 
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where the overall normalization of 2% is chosen to make contact with AMY conventions. 
Each component of / satisfies the same Schrodinger equation (2.8) that G does, so that 



idtf(B,t) = Hf{B,t) 



with initial condition 



f(B,0) = -2 i V B 5^(B). 
Now define the time-integrated amplitude 



f(B) 



dtf(B,t). 



(A8) 
(A9) 

(A10) 



Integrating both sides of (A8) over time (and noting that f(B,t) decays with time because 
of the — iT 3 piece of H), 

-2V B 5 (2 \B)=Hf(B). (All) 
The rate (A6) can be written in terms of f(B) as 

cT'brcm (XX P s —+g(x^ 



00 



duo [x(l-x)Ef 



Re 



(2t)- 1 V B -/(B) 



B=0 



(A12) 



Now Fourier transform from B to p B . Using the form (2.7) of H, the equation for / becomes 
- 2ip B = 5E(p B ) f(p B ) - i J d 2 q ± §^\\C A [f(p B ) - f(p B + q ± j 

+ (C s - \C A ) \f(p B ) - f(p B + xq ± )] + \cAf{p B ) - f{p B + (1 - x)q ± ) 



(A13) 



Instead of p B , AMY uses the variable h = p B P. In the case of bremsstrahlung, they define 
the momenta of the splitting particles as p' = P, k = xP, and p — (1 — x)P. If one defines 



F(h) = Pf(h/P), 



(A14) 



then (A13) becomes 



-2ih = 5EF(h)-i I d 2 q ± dTd J 



d 2 q A 



C A 



F(h)-F(h + p'q ± ) 



+ (C s - \C A ) F(h) - F(h + kq ± ) 



F{h)-F(h + pq±) 



(A15) 



This is equation (6.7) of Ref. [12] if one changes integration variable from q± to — q± in some 
of the terms and recognizes that 



d 2 q± (2tt) 



2 2 

9 A( \- 9 



/+oo 
-oo 



d £{A-{Q)[A-{Q)\ //q ^- 



(A16) 



With the same notation, the rate (A12) becomes 



^-f^brem ^-f^brem 

CO — = x ■ 



ax P s ^ g (x) 



duo 



dx 



d 2 p B 

4[x(l - x)E} 2 J (2tt) 2 



Re[2p B ■ f(p B )\ 



Ax 2 (l-x) 2 E & J (2tt) 



j-0- 2 Re[2h-F(h)}. (A17) 
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This formula can be extracted from the rates per unit volume presented for kinetic theory 
in AMY Ref. [13], for example, with 

W£; & - = {2n ^J^% E s J <gl ^ • F( h )]. (A18) 

More simply, if final-state factors of [1 ± f{xE)\ [1 ± /((l — x)E\ are included, it corresponds 
to Eq. (5) of Jeon and Moore [15] or Eqs. (1.1) and (4.1-2) of Ref. [19]. Jeon and Moore 
use the symbol dT/dt to denote rate rather than T. 

Readers comparing to AMY should beware that AMY uses the symbol T to indicate the 
rate per unit volume, integrating what I call T^, rem over the initial particle's momentum with 
a factor of its distribution function / and including final state factors. 17 



Wiedemann 



Finally, I will translate to the notation of Wiedemann and collaborators [22, 24-26] as 
presented in Salgado and Wiedemann [25]. They specialize to the x € 1 limit of soft 
bremsstrahlung gluons, but they study more properties of the process, such as the angle 
between the emitted gluons and the high-energy parton, and what happens when the gluon 
momentum is so small that the approximation k±_ <C k is no longer valid. The basic result, 
Eq. (2.1) of Ref. [25], is 



dl 



s^R 



du (27r) 2 c<r 



2 Re 



dyi 



Co 



III 



dyi du dk±e 



x e -(l/2)/~den(0<r(«) 



d_ d_ 
dy du 



x 



Jy=0 



u=r(y L ) 



Vr exp 



.<// 



m 



s 2 V ioo 



(A19) 



The limit k±_ < is used to restrict attention to gluon bremsstrahlung in a finite opening 
angle G with x — sinO. In this paper, I have put no such restriction, and I have assumed 
k sufficiently large that k± k dominates. This corresponds to replacing the upper limit 
Xu on the kj_ integration by infinity. That integral then generates a factor of S^(u), which 
makes the u integration trivial. Using the fact that their definition of a(u) has a(0) = 0, 
one then obtains 



u ■ 



dl a s C R 



duj 



2 Re dyi / dy, 



x 



d_ d_ r= 

dy du J v=r , 



u=r(y L ) 



Vr exp 



in 



s 2 



n(0 a{r) 



iuj 



(A20) 



u=y=0 



17 A pernicious factor of 2 that arises when comparing to AMY expressions is that they sum formulas for 
splitting of particle types a — > be over the types b and c. For bremsstrahlung from a quark, this gives rise 
to a factor of 2 because both q — > qg and the identical q — > gq are summed over. For g — > gg there is no 
such factor of 2, accounting for the relative factor of 1/2 one needs to include when integrating over the 
final momentum fractions of two identical particles. 
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J_1 ■ 

this paper 


Salgado & Wiedemann [25] 
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Vi 






m 






y 




B 2 


u 




B 


r 




q 


q/C A 




M 


UJ 










* 


uj L 






* 


r 3 (B) 


no(r) 
2 




f 2 (B) 


n cr(r) 
2Ca 





TABLE II: Translation between notation of this paper and Salgado and Wiedemann [25], which 
studies the i«l limit. 

This is the small x approximation to (2.6) of this paper, with the notational translations 
shown in table II, my convention £ — 0, the bremsstrahlung gluon mass ignored, and the 
Green function expressed as a path integral. 
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